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Abstract 

This paper gives a simple and original presentation of various coronagraphs inherited from the 
Lyot coronagraph. We first present the Lyot and Roddier phase mask coronagraphs and study their 
properties as a function of the focal mask size. We show that the Roddier phase mask can be used 
to produce an apodization for the star. Optimal coronagraphy can be obtained from two main 
approaches, using prolate spheroidal pupil apodization and a finite-size focal mask, or using a clear 
aperture and an infinite mask of variable transmission. 

1 Introduction 

Direct imaging of faint companions or planets around a bright star is a very difficult task, where the 
contrast ratio and the angular separation are the observable parameters. The problem consists of 
detecting a faint source over a bright and noisy background, mainly due to the diffracted stellar light. 
High contrast ratios and small angular separations correspond to the most difficult case. Typically, 
for giant exoplanets, contrast ratios of about 10~ 7 are expected in the near infrared (J;H;K bands), 
based on models for relatively young objects of about 100 Myr [71 SI E] . According to these models, 
older objects would be an order of magnitude fainter. Terrestrial planets are much fainter than giant 
planets, about 3 to 4 orders magnitudes fainter depending on the wavelength range. 

The aim of this article is to give a short review of basic concepts and techniques used in focal plane 
mask coronagraphy for exoplanets detection. The paper will focus on derivation of a general formalism 
which allows to gain deeper insight in the behaviour of the principal coronagraphic techniques. It will 
not consider important problems related for example to the effects of adaptive optics residual speckles, 
slow-varying speckles caused by mechanical or thermal deformations, telescope central obstruction, 
chromaticity. 

2 Coronagraph general formalism 

A common formalism can be used for Lyot and Roddier & Roddier coronagraphs [15, 19J in their 
classical or apodized version j3j |23l [22] , four quadrants coronagraph (4QC) [20 , band-limited mask 
coronagraph p3] and shaped pupil coronagraph [12] . 

*R. Soummer is supported by a AMNH Kalbfleisch Research Fellowship. 
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Figure 1: Schematic illustration of Lyot, Roddier & Roddier, 4QC and band-limited mask coronagraph. 



We will follow the notations of [3]. The successive planes of the coronagraph are denoted by A, B, 
C and D. A is the entrance aperture, B denotes the focal plane with the image mask, C is the image 
of the aperture with the Lyot stop and D is the image in the focal plane after the coronagraph. The 
general setup is illustrated in figure [TJ The position vector in each plane of the coronagraph will be 
noted in bold. 

• We denote by p(x) the telescope pupil and a(x) its possible transmission if a pupil apodization 
is used. We assume that x in A is in units of wavelength. 

• The mask transmission is t{x). We will consider without loss of generality that t{x) = 1 — mix). 
The units for the coordinates in B and D are angles on the sky. 

• The Lyot stop, denoted as p s (x) will clearly acts as a low pass filter and consequently its size 
will be chosen at most equal to the size of the entrance pupil so that p s (x)p(x) = p s (x). 

We denote as V, S and A4 the areas defined by the pupil, the stop and if pertinent the mask. We will 
make the usual approximations of paraxial optics. Moreover we neglect the quadratic phase terms 
associated with the propagation of the waves or assume that the optical layout is properly designed 
to cancel it [1. The constant propagation terms between the successive planes will be omitted. In 
this case the coronagraph can be easily described using classical Fourier optics: a Fourier transform 
exists between each of the four planes. 

The wavefront complex amplitude just before the pupil plane is ^(x). The expression of the 
complex amplitude in the successive planes A, B, C and D are: 



^a(x) = ^(x)p{x)a{x) 

^b(x) = t(x)i>A(x) = ^a(x) - m(x)$A(x) 

^c(x) = ^>b{x)Ps{x) = (*a(£) - ^a(x) * rh(x))p s (x) 

^ D (x) = ^c(x) = ($ A (x) ~ $A{x)m(x)) *p s (x) 



(1) 
(2) 
(3) 
(4) 



where g is the Fourier transform of g and * denotes convolution. 

The effect of the coronagraph clearly appears in equation ([3]). The first term is the direct wave 
diaphragmed by the Lyot stop. The second term corresponds to the wave diffracted by the mask 
for which the light diffracted outside the aperture in C has been removed. A coronagraph correctly 
designed for exoplanets imaging can operate one of these two techniques: 

1. cancel the on-axis star without altering an off-axis source: for the star the two terms in Q must 
interfere destructively whereas for the planet the second term of ([3]) must cancel. In this case 
the response of the coronograph is no longer invariant by translation. 

2. concentrate the on-axis star light reducing the off-axis diffracted light. This approach corre- 
sponds to the apodization techniques [1]. 
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These two techniques can be achieved by a proper choice of the apodization a{x) and the transmission 
t{x). This article will focus on the first solution. For a general overview of pure apodization techniques 
see [T] and included references. 

Section 2 will present the main solutions to this problem that have been derived in the litterature. 
First the "historical" Lyot coronagraph and the Roddier coronagraph are presented. Whereas the first 
attempts to cancel the star light, it appears that the second one can also operate by an apodization 
of the star light. Whereas these coronagraph can provide acceptable results when the dynamic is not 
too high their performances appear to be insufficient for the detection of earth-like planets. Section 
3 addresses the problem of optimal coronagraphy which try to maximize a criterion quantifying the 
star light rejection. In this context a complete star light rejection can be obtained using an infinite 
size mask (4QC and band-limited mask) or a finite size ir phased mask (Roddier coronagraph) with a 
properly apodized entrance pupil. In the case of a Lyot coronagraph, an apodized aperture allows to 
maximize the star light rejection for a given mask size. Through all this presentation we will try to be, 
when possible, as general as possible with respect to the geometry of the system without specializing 
for example on a specific pupil shape. 

3 Lyot and Roddier coronagraphs 

The first solution was proposed by Lyot |15| and consists in canceling the major contribution of the 
star energy located in the telescope point spread function (PSF) inside a disk of radius r m . This is 
simply achieved setting for t{x) an opaque mask at the center of the image plane B. 

The results will be presented for a circular aperture of radius R and without central obstruc- 
tion. The use of polar coordinates will be preferred. Under this asumption p{x) = H{r/R), m{x) = 
n(r/ r m) and p s (x) = II(r/r s ) where Il(r) = 1 for r 6 [0,1) and if r > 1 and for example: 
p{x) = Rr^ 1 Ji(2-7ri?r). Note that as long as ^a{x), a(x) and m(x) are radial functions, the complex 
amplitudes at the different stages of the coronagraph will exhibit the same symmetry. 

Figure [2] illustrates the response of the Lyot coronagraph to an on-axis point source (ty(x) = 1) 
for different values of the mask radius r m . These values of r m are located on ^^(r) in figure pi 

We seek to obtain the best subtraction of the two wavefronts ^ a{t) and ^A( r ) * rh{r) inside de 
Lyot stop. Figure [2] cleary shows that this result is achieved increasing the mask size. In fact, as 
r m increases, rh{x) will be more "concentrated" around the origin. As long as rh{x) always verifies 
f m{x)dx = m(0) = 1 we have m{x) —* 5{x) as r m — > +00, which is the neutral element for the 
convolution product. On the contrary as r m decreases m{r) widens around m(0) = vrr^. At the limit 
r m —>■ 0, m(x) — ► 7rr^ and ^a(x) * m(x) tends to the constant ^R 2 ^, which does not substract to 
^a(x) = 1. It is interesting to note that this problem is formally equivalent to the problem of digital 
low-pass filter design with finite impulse response using the windowing method where the infinite 
impulse response of the ideal filter is truncated using a window, see for example |17j. 

The results of figure [2] suggest the following comments: 

• The performance of the coronagraph is not strictly linear with the mask size. The choice of 
mask b is as an evidence more appropriate than mask c. This is confirmed by the general shape 
of the second plot in figure [3] which gives the integrated intensity in C for a Lyot stop of radius 
R. 

• A large amount of the residual star energy is located at the edges of the pupil. Consequently a 
moderate reduction of the diameter of the Lyot stop radius r s will provide a significant gain in 
starlight rejection with a reasonable loss of transmission. The effect of a Lyot stop reduction is 
visible in the last two lines of the plot. The value of r s for the last line is indicated in the plots 
of the intensity in C by a dashed line. 

• The off-axis transmission of the coronagraph is of course a crucial point for planet detection. As 
long as the angular distance from the optical axis is sufficient to guaranty that the response of 
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the planet in B will not interfere with the mask, we can consider that the response of the planet 
will be the shifted PSF of a pupil with radius r s . In order to quantify the loss of transmission 
for an off-axis source, this PSF has been added in figure [2} Finally, note that a solution to the 
analytical computation of ^c{ r i Q)i *&D( r , 6) of a point source which can be close to the optical 
axis has been addressed in [8] . 



Lyot mask a 



Lyot mask b 



Lyot mask c. 

















\ 


/-o 


.5 





5 



^0.5 











4 










.5 


0.5^ 





0.25 






0.25 






0.25 




— 






I 




J 




li / 


LJ 


-0.5 




0.5 


-0.5 




0.5 


-0 


5 


0. 




2 14 6 8 10 






Figure 2: Illustration of Lyot coronagraphy for circular aperture (radial cuts) and an on-axis point 
source. The first raw represents the two complex amplitudes that should interfere destructively in 
C. The second row gives the intensity in C and the third and fourth raw the intensity in D. In the 
two last rows, the continuous line represents the intensity in D for a coronagraphed point source and 
the dashed line the PSF associated to a pupil of radius r s . The vertical dashed line in the last row 
indicates the radius of the mask. 

The Lyot technique was improved by Roddier [19] replacing the opaque mask by a tt phase mask: 
m(x) = 2U(r/r m ). Figure [4] show that for very small values of r m , the Roddier coronagraph tends 
to behave like the Lyot coronagraph, trying to cancel out the light of the star. This behavior totally 
differs when r m increases: for mask size a or b the effect of this coronagraph is to produce in the plane 
C an apodized version of the wavefront originated by an on-axis source. Note that a coronagraphic 
technique relying on a similar principle was proposed in [16 a using a ir/2 phase mask and a defocus. 
The effect of this apodization is visible in D: the main lobe of the response is broadened whereas the 
side lobes decrease. Moreover, it is worthy to note that this apodization as a maximum value for x = 
that can be higher than 1. 

These results suggest the following comments: 

• the choice of r m remains a crucial point. However, according to the preceding remark the 
integrated energy in C ploted in figure [3] is no more a valid criterion for selection of an optimal 
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Figure 3: The left plot locates the different masks used in the simulation on the response of an on- axis 
point source. The right plot gives the integrated intensity in C for r s = R as a function of the mask 
size. 

value of r m . 

• As the second row of figure [3] shows, a reduction of the Lyot stop is not crucial as long as the 
Roddier coronagraph operates in its " apodization mode" . 

• The effect of a Roddier coronagraph on a planet relies on the same principle as the Lyot coro- 
nagraph. Consequently the apodization effect mentioned above will not operate for an off-axis 
source. Moreover, the fact that for this latter coronagraph the mask will be smaller and a 
reduced stop is not necessary will increase its performances in terms of resolution and planet 
flux. 

According to the previous remarks, a Roddier coronagraph with an extended mask can be used as an 
apodizer for the wavefront complex amplitude. The star rejection of such a system being insufficient 
for planet detection such a device should be considered as the first stage of a "classical" coronagraph. 



4 Optimal coronagraphy 

According to equation ([3| the ideal coronagraph equation is: 

#a(x) = ^a(x) * m(x), Vf e S (5) 

where ^a(x) = a{x)p{x). Following Lyot and Roddier coronagraphs, a first solution is to set m{r) = 
eIL(r/r m ) where e = 1 for Lyot and e = 2 for Roddier and to find, if it exists, the optimal apodization 
a(x) solution of the integral equation (|5]). This solution was first proposed by [5] and [11]. The other 
solution is to suppress the constraint on the mask and to find the function m{x) solution of equation 
^ when a(x) = 1. This approach includes the band-limited masks proposed by [2] and the four 
quadrants coronagraph |20j . 



4.1 Finite size masks: apodized Lyot and Roddier coronagraphs 

As long as m{x) has a bounded symmetric support Ai, we can write using the notations of Apendix 
B: m(x) = €Km(x). Consequently, equation (13) implies that if a(x) is proportional to a spheroidal 
prolate function associated to Ai and V, the residual in C when p s {x) = p{x) is obtained replacing in 
^ the second term by eA^^(x): 

<S c (x) = (1 - eA)^(z) (6) 

which states that the residual wavefront in C originated by an on-axis source is proportional to the 
apodized entrance pupil. If we assume that the angular distance of the planet is sufficient, its response 
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Figure 4: Illustration of Roddier coronagraphy for circular aperture (radial cuts) and an on-axis point 
source. The first raw represents the two complex amplitudes that should interfere destructively in C. 
The second row gives the intensity in C and the third raw the intensity in D. In the last row, the 
continuous line represents the intensity in D for a coronagraphed point source and the dashed line the 
PSF associated to a pupil of radius r s . 

in C is obtained setting e = in (p]). The tradeoff between the reduction of the planet flux by the 
apodization and the reduction of the star flux by the coronagraph can be evaluated by the quantity 
rj defined as the ratio between the total planet intensity in D and the total star intensity outside the 
mask in D. 

From equation (j(3| the star intensity in D is simply (1 — eA) 2 |^^(af)| 2 . The total energy outside 
the mask is: 

(l-eA) 2 {[ \V A (x)\ 2 dx- I \^> A (x)\ 2 dx) =(1-6A) 2 ( f \* A (x)\ 2 dx - A I \* A (x)\ 2 dx) (7) 
Jm. 2 Jm Jul 2 Jm. 2 

where the last equality comes from the fact that eA is the ratio between the energy encircled in the 
bounded frequency domain A4 and the total energy of ^ A (x). As a consequence we have: 

r ? = (l-eA) 2 (l-A) (8) 

The objective will be to minimize r\ for a given coronagraph type (e = 1 or 2). This is achieved by a 
modification of the size of M. which will affect A and the pupil apodization shape a{x). 

Any eigenvalue of the integral equation minimizing ([8]) leads to a valid solution. However we will 
only consider in the sequel the largest eigenvalue, i.e. corresponding to the maximum encircled energy 
behind the focal plane mask. This choice will in fact "generally" lead to a positive prolate function 
which will be normalized by its maximum value on the pupil in order to achieve an apodization with 
maximum transmission. Finally, it is worthy to note that this development does not require a specific 



shape of the pupil and the mask. However, further analytical derivations based on equations (14) 



and (15) require that the mask is a scaled version of the pupil: in this case if the pupil has a central 



obstruction m(x) has a hole. This point will not be considered herein. 
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• We consider first the case of an apodized Roddier coronagraph, rj = (1 — 2A) 2 (1 — A). The results 
of Appendix A prove that 1/2 is a valid eigenvalue and consequently a mask size can be chosen 
to achieve rj = 0, i.e. a total extinction of the star. 

For example, in the case of a circular pupil without obstruction, figure [8] shows that the coefficient 
c = r m R ~ 0.25 corresponds to A = 1/2. This value defines the apodization shape and fixes 
the mask size to r m k, 0.25/ R. Figure gives the corresponding apodization. This figure also 
illustrates the loss in transmission for an off axis planet comparing the prolate apodized pupil 
PSF with the PSF of the unapodized pupil. 

• We consider now a prolate apodized Lyot coronagraph (PALC), rj = (1 — A) 3 . The eigenvalues 
being upper-bounded by 1, a prolate apodized Lyot coronagraph (PALC) cannot achieve total 
extinction of an on-axis point source. The trivial solution A = 1 would correspond to an infinite 
size opaque mask. However, approximate solutions can be obtained for eigenvalues A close to 1 
and finite mask size. Taking advantage of the rapid saturation of the eigenvalue curve we can 
choose a mask size corresponding to an eigenvalue close to 1. 

Figure [6] illustrates the response of the PALC to an on-axis point source for different values of 
the mask radius r m given in figure [3j To each value of r m corresponds a coefficient c and and 
optimal apodization function given in the first row. The second row gives the residue in C and 
the last row the response of the star and the PSF associated to the apodized telescope pupil. 
Note that contrarily to the unapodized Roddier coronagraph operating in its "apodized regime" 
(see column 2 and 3 of figure [4]) the residue in C is here both apodized and attenuated. An 
analog comment can be done for the apodization alone coronagraphs. 

The property that the Lyot stop wave amplitude is proportional to the entrance apodized pupil 
amplitude, creates the possibility for multiple stage coronagraphs described in [2]. A multiple 
stage PALC only requires a single apodizer in the entrance pupil. The Lyot stop plane is naturally 
apodized and can be used as the entrance pupil of a second coronagraphic stage without further 
loss of throughput due to apodization. In this case the residual amplitude in the second Lyot 
stop plane is (1 — A) 2 ^a(x) and the parameter 77 becomes (1 — A) 5 . 




-0.5 0.5 -1 1 



Figure 5: Illustration of Roddier apodized coronagraphy for circular aperture (radial cuts). The left 
plot shows the optimally prolate apodized pupil. The right plot shows the PSF corresponding to the 
prolate apodized pupil and the PSF of the unapodized pupil. 



4.2 Infinite size masks band limited and four quadrants coronagraphs 

A simple solution to define a mask m(x) that verifies ^ is to chose m(x) such that rh(x) is non zero 
only on a small bounded region B included inside the pupil. In this case if we define S such as Vxo £ <5 
the support of rh{xQ — x) is strictly included in the pupil V: 

Va?o € 5, ^a(^g) * m(xo) = / rh(xQ — x)dx = m(0) (9) 

Jv 
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Figure 6: Illustration of prolate apodized coronagraphy for circular aperture (radial cuts) and an 
on-axis point source. The first raw represents the two complex amplitudes that should interfere 
destructively in C. The second row gives the intensity in C and the third and fourth raw the intensity 
in D. In the last row, the continuous line represents the intensity in D for a coronagraphed point 
source and the dashed line the PSF associated to the apodized telescope pupil. 

Consequently, if we impose m(0) = 1 we will clearly verify equation ^ on S C V. Figure [7] illustrates 
the definition of the various sets S, V and B. A simple solution to block the residue in V\S is to 
use as a Lyot stop the indicator function of S, p s (x) = ts(%) (Ps(%) = 1 if x € 5 and elsewhere). 
Consequently a properly defined band limited, and necessary infinite size mask m{x) must verify: 
m(0) = 1, < m(x) < 1 and m{x) is bandlimited on B C V. This approach was proposed by [15] . 

A fundamental result for understanding of band limited coronagraphs states that the complex 
amplitude in D associated to an off-axis point source at an angle x\ is the response of the Lyot stop 
attenuated by the transmission at x\\ 

^ D (x) =t(x 1 )p s (x-x 1 ) (10) 
This result is demonstrated in [13]. This result suggests the following comments: 



1 . The loss of throughput of a bandlimited coronagraph is principally due to the undersizing of the 
Lyot stop. Consequently B should be as small as possible. However the general consequences of 
a reduction of B are a widening of the central lobe of mix) and an increase of the ripples in the 



tails of m(x). From equation (10) the first effect implies an increase of a blind zone arround the 



star and the second a reduction of detectability for given locations of a planet. 
2. The behavior of the coronagraph when the star is slightly off- axis is perfectly described by 



equation (10): in order to reduce the effect of a misalignmemt m(x) must be as "flat" as possible 
in 0. This effect is quantified in [13 j by the degree of the first term in the serie expansion of 
t(x). Note that the flatness of m{x) around x = is of course deeply related to the constraints 
previously mentioned. Finally, a technique to construct a mask of given order adding and 
multiplying simple band-limited masks is proposed in [13] . 
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Figure 7: Left figure: relation between the pupil (V), support of the mask Fourier transform (B) and 
the Lyot stop S in a bandlimited coronagraph. Right figure: illustration of three different masks. 
Dashed lines: second and fourth order Bessel cardinal masks. Continuous line: prolate spheroidal 
function. 

To conclude, figure [7] compares the behavior of three different radial transmissions in the case of a 
circular aperture. To facilitate this comparison the three masks have been normalized to the same B in 
order to guaranty that the Lyot stop will be the same: for these simulations the flux reduction equals 
4/9. The first transmission is a Bessel cardinal function and the second one the square of a Bessel 
cardinal function. These functions have been properly scaled and shifted in amplitude in order that 
t(r) occupies all the interval [0, 1]. The properties detailed in the item (i) of the discussion above can 
also be expressed in term of concentration of the energy of m{x) around the origin. In this case, the 
additional band-limited constraints leads naturally to a prolate function. Such a solution is illustrated 
in figure [7} In this case the coefficient c has been chosen empirically equal to 6.8 in order to achieve 
a good compromise between the the ripples and the width of the main lobe. This result shows that 
the behavior of such a prolate mask and a squared Bessel cardinal mask are similar. Finally, it is 
important to emphasize that the masks illustrated in figure [7] correctly behave with respect to the 
requirement (i). On the contrary they do not fulfil (ii) (the order of the squared Bessel cardinal and 
the prolate is only 4) and will have a poor behavior in the case of a resolved star or a misalignment. 

To conclude, it is important to mention in this section devoted to optimal coronagraphy with 
infinite size masks the 4QC. This coronagraph will be developed in a separate paper by D. Rouan in 
this volume. It relies on the fact that for a circular aperture and a transmission t{x) = sign(rr)sign(y), 
the complex amplitude in C is identically zero inside the pupill for a point source on the axis. This 
beautiful result relies on a nontrivial property of the Fourier transform (see the paper of D. Rouan and 
included references for a proof) . A coronagraphic technique based on a similar principle is developped 
in [9]. 
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A Prolate spheroidals functions 

This section presents some facts about prolate spheroidals functions. For a detailed presentation refer 
to the seminal papers [21] and references therein. Application of prolate spheroidals functions to optics 
can be found in [10J. 

In order to gain deeper insight in their principal properties using simple mathematics the cartesian 
coordinates have been prefered. Derivations using specific coordinates can be found in the references. 
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We consider a real valued square-integrable function / of d variables having a bounded support V (in 
our case f(x) = a(x)p(x) and V is the telescope pupil): 



f(x) = / f(y)e-^dy 
Jv 



(11) 



The ratio between the energy encircled in the bounded frequency domain A4 (here the Lyot mask) 
and the total energy of / is obtained integrating the squared modulus of (11 ) on A4: 



K M (x - y)f(y)f(x)dydx/ / f(y) 2 dy 
v Jv 



(12) 



where K_m(x) is the inverse Fourier transform of the indicator function of M, H_m(x). Standard 
results on functional analysis, see for example [18], prove that the maximum of this ratio is the largest 
eigenvalue A of the integral equation: 



v 



K M {x - y)ip(y)dy = Aip(x), x G V 



(13) 



or equivalently: (tp(x)l.p(x)) * Km(x) = A-0(x), x G V. 

The prolate spheroidal functions are defined as the solution of this integral equation. The extrema 



of (12 ) is reached when f(x) is the corresponding eigenfunction. The kernel being positive defined the 



eigenvalues of (13) are positive and from (12) upper bounded by 1. The eigenfunctions of (13) are 
orthogonal and complete on L2(V) and orthogonal on M. d when extended outside V using (13). 



Considerable simplifications occur when V and A4 are both scaled versions of a "normalized" 
domain U: 3R G R d , 3f m G M d such that x G V 44> x./R eU, x G M O x./f m G U, where a./b (resp. 
a.b) defines the vector having a(k)/b(k) (resp. a(k)b(k)) as components. A simple change of variable 

in the definition of K_m(%) shows that K_m(x) = (jYl=\ r m {k)^ Ku(r m .x). Substitution of this result 

(j)(R.x) where q 



in (13) shows that ip(x) 



is the solution of the normalized problem: 



J[r m (k)R(k)J J^K u (? m .R-(x 



y))(j>{y)dy = A0(x), x£U 



(14) 



Consider for example the case where d = 2 and U is a disk. Consequently V and M. are ellipses. 



Equation ( 14 ) shows that the prolate functions associated to an ellipse are obtained by a scaling 



of the prolate function associate to a disk. 



• If Vfc, R{k) = R and r m {k) = r m , equation (14) shows that the prolate spheroidal functions 
for a given U form essentially a one-parameter family of functions depending of the parameter 
c = r m R. 

Additional simplifications in this last case are obtained when the region IA is symmetric: x £U 44> 



—x G U. In this case the solutions of ( 14 ) are real, either even with a real eigenvalue or odd with a pure 



imaginary eigenvalue. Moreover, finding the solutions of (14) is equivalent to finding the solutions of: 



am [xj 



-i2ircx t 



y (j)(y)dy, iGW, where A 



\a\ 2 c d 



(15) 



u 



This equation states that the Fourier transform of 4>{x)\i{x) is proportional to 4>(x/c). This result 
can be easily verified taking the complex conjugate of (|15[), substituting in the integral q^(y) by the 



right side of (15) and identifying with (14) 



Circular prolate functions correspond to the case d = 2 where U is the unit disk. They decompose, 
up to a normalisation factor as 4>N,n(r,9) = i?jv,n(^)[cos, sin](A r 0). Consequently (15) shows that they 
can be defined by their invariance to a finite Hankel transform of order N, e.g.: 



l 

2vr / Ro.o(r)Jo(2irrp)rdr 



A 



R 



0,0 



P 



(16) 



10 



Figure [8] represents the radial circular prolate function </>o,o( r ) f° r three different values of c. 

A crucial point is the numerical computation of prolate functions. For simple geometries, the func- 
tions can be computed using rapidly converging series as the one derived for a circular aperture in [21] ■ 



For general geometries, prolate spheroidal functions can be computed directly solving equation (13) 



using one of the numerous iterative algorithm that can be found in the litterature for the computation 
of the eigenelements of linear integral operators, e.g. the algorithm used in |llj . 




Figure 8: Circular prolate functions. Left plot gives the largest eigenvalue A as a function of 4c. The 
right plot gives the prolate functions for three values of c. 
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